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ABSTRACT 

The compact, nonthermal radio source Sgr A* at the Galactic Center 
appears to be coincident with a ~ 2.6 x 10 6 M & point-like object. Its energy 
source may be the release of gravitational energy as gas from the interstellar 
medium descends into its deep potential well. However, simple attempts at 
calculating the radiative spectrum and flux based on this picture have come 
tantalizingly close to the observations, yet have had difficulty in accounting for 
the unusually low efficiency in this source. Regardless of whether the radiating 
particles in the accretion flow are thermal or nonthermal, there now appear to 
be two principal reasons for this low conversion rate of dissipated energy into 
radiation: (1) the plasma separates into two temperatures, with the protons 
attaining a significantly higher temperature than that of the radiating electrons, 
and (2) the magnetic field B is sub-equipartition, which reduces the magnetic 
bremsstrahlung emissivity, and therefore the overall power of Sgr A*. In this 
paper, we investigate the latter with a considerable improvement over what has 
been attempted before. In particular, rather than calculating B based on some 
presumed model (e.g., equipartition with the thermal energy of the gas), we 
instead infer its distribution with radius empirically with the requirement that 
the resulting spectrum matches the observations. Our assumed ansatz for B(r) 
is motivated in part by earlier calculations of the expected magnetic dissipation 
rate due to reconnection in a compressed flow. We find reasonable agreement 
with the observed spectrum of Sgr A* as long as its distribution consists of 3 
primary components: an outer equipartition field, a roughly constant field at 
intermediate radii (~ 10 3 Schwarzschild radii), and an inner dynamo (more or 
less within the last stable orbit for a non-rotating black hole) which increases B 
to about 100 Gauss. The latter component accounts very well for the observed 
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sub-millimiter hump in this source. 

Subject headings: accretion — black hole physics — hydrodynamics — Galaxy: 
center — magnetic fields — magnetohydrodynamics — plasmas — turbulence 
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1. Introduction 

The Galactic center (GC) has long been suspected of harboring a central mass 
concentration, which appears to be coincident with the unique, nonthermal radio source, 
Sgr A*. Haller et al. (1996) used the velocity dispersions of stars at ^ 0.1 pc from Sgr 
A* to derive a compact mass of ~ 2 x 10 6 M Q . This is consistent with the value of 
~ 2.5 — 3.2 x 10 6 M Q derived more recently by Genzel et al. (1996), using the radial 
velocities and velocity dispersions of ~ 25 early-type stars and of ~ 200 red giants and 
supergiants within the central 2 pc. A third technique for tracing the central gravitational 
potential is based on the acquisition of proper motions for the ~ 50 — 100 brightest stars 
within the radial range ~ 0.004 - 0.4 pc (Eckart & Genzel 1996; Ghez et al. 1998). These 
stellar motions seem to require a central dark mass of (2.6 ± 0.2) x 10 6 M , in good 
agreement with earlier ionized gas kinematics and the velocity dispersion measurements. 

Of course, showing that the GC must contain a centralized mass concentration does 
not necessarily imply that this dark matter is in the form of a compact object with a few 
million solar masses. It does not even imply that the unusual radio source Sgr A* must 
be associated with it. VLBA images of Sgr A* with milliarcsecond resolution (Bower & 
Backer 1998) show that at A7 mm, its size is 0.76 ± 0.04 mas, or roughly 6.2 x 10 13 cm, 
much smaller than ~ 0.01 pc, the present limiting region within which the 2 — 3 x 10 6 M Q 
are contained. So the dark matter may be distributed, perhaps in the form of white dwarfs, 
neutron stars, or ~ 10 M Q black holes (e.g., Haller et al. 1996), though the latest stellar 
kinematic results appear to rule out the first two possible constituents (Genzel et al. 1996). 

Whatever the composition of a distributed mass concentration is, one is left with 
the task of accounting for the nature of Sgr A* itself. It is likely that many of Sgr A*'s 
characteristics are associated with the liberation of gravitational energy as gas from the 
ambient medium falls into a central potential well (Melia 1994; see also Ozernoy 1989 for an 
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alternative conclusion regarding wind accretion). There is ample observational evidence in 
this region for the existence of rather strong winds in and around Sgr A* itself (from which 
the latter is accreting), e.g., the cluster of mass-losing, blue, luminous stars comprising 
the IRS 16 assemblage located within several arcseconds from the nucleus. Measurements 
of high outflow velocities associated with IR sources in Sgr A West (Krabbe et al. 1991) 
and in IRS 16 (Geballe et al. 1991), the H 2 emission in the circumnuclear disk (CND) 
from molecular gas being shocked by a nuclear mass outflow (Genzel et al. 1996; but see 
Jackson et al. 1993 for the potential importance of UV photodissociation in promoting this 
H 2 emission), broad Bra, Br7 and He I emission lines from the vicinity of IRS 16 (Hall et 
al. 1982; Allen et al. 1990; Geballe et al. 1991), and radio continuum observations of IRS 7 
(Yusef-Zadeh & Melia 1991), provide clear evidence of a hypersonic wind, with a velocity 
v w ~ 500 — 1000 km s -1 , a number density n w ~ 10 3 ~ 4 cm -3 , and a total mass loss rate 
M w ~ x 10( _3 )~( -4 ) M , pervading the inner parsec of the Galaxy. Even so, the observations 
do not yet provide sufficient information for us to identify the physics of accretion when the 
infalling gas penetrates to within about 10 3 or 10 4 Schwarzschild radii of the central object. 



1.1. Behavior of the Accreting Gas at Small Radii 

Three-dimensional hydrodynamic simulations (Coker & Melia 1997) indicate that 
the accreted specific angular momentum A (in units of cr s , where r s = 2GM/c 2 is the 
Schwarzschild radius in terms of the black hole mass M) can vary by 50% over £ 200 
years with an average equilibrium value for A of 40 ± 10. Thus, even with a possibly 
large amount of angular momentum present in the wind, relatively little specific angular 
momentum is accreted. This is understandable since clumps of gas with a high specific 
angular momentum do not penetrate to within 1 Ra, where 

R A = 2GM/v w 2 (1) 
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is the capture radius defined in terms of the wind velocity v w at infinity. The variability in 
the sign of the components of A suggests that if an accretion disk forms at all, it dissolves, 
and reforms (perhaps) with an opposite sense of spin on a time scale of ~ 100 years. 

The captured gas is highly ionized and magnetized, so it radiates via bremsstrahlung, 
cyclo-synchrotron and inverse Compton processes. However, for purely spherical accretion, 
the efficiency of converting gravitational energy into radiation is quite small (as little as 
10~ 4 in some cases), so most of the dissipated energy is carried inwards (Shapiro 1973; 
Ipser & Price 1977; Melia 1992). In fact, if the magnetic field is a negligible fraction of its 
equipartition value (see below), Sgr A* would be undetectable at any frequency, except 
perhaps at soft X-ray energies. But as the plasma continues to compress and fall toward 
smaller radii, one or more additional things can happen, each of which corresponds to a 
different theoretical assumption, and therefore a potentially different interpretation. 

The questions one may ask include the following: (1) Does the flow carry a large 
specific angular momentum so that it forms a disk with lots of additional dissipation? (2) 
Does the flow produce a radiatively dominant non-thermal particle distribution at small 
radii (e.g., from shock acceleration), or does thermal emission continue to dominate the 
spectrum? (3) Does the flow lead to an expulsion of plasma at small radii that forms a 
non-thermal jet, which itself may then dominate the spectrum? These, either individually 
or in combination, have led to a variance of assumptions about the nature of the inflowing 
gas that then form the basis for the development of different interpretations. 

Observationally, a key issue is why the infalling gas maintains a low radiative 
efficiency. Beckert & Duschl (1997) suggest that shocks in the accreting plasma produce 
a power-law electron distribution, which is truncated by strong cooling. This forms a 
"quasi" mono-energetic distribution. The overall emission, which is strictly non-thermal, is 
suppressed by constraining the number density of relativistic particles and the intensity of 
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the magnetic field (at about 5 to 10 Gauss). 

Falcke, Mannheim and Biermann (1993) and Falcke & Biermann (1999), on the other 
hand, assume that the infalling plasma eventually produces a jet of power-law electrons 
whose number density varies with radius in the expulsion. The overall emission, which is 
a sum of non-thermal components, is also suppressed by constraining the particle number 
density and hence the equipartition magnetic field, both of which are assumed to be scaled 
by a low luminosity disk. 

In the picture developed by Narayan, et al. (1998), the infalling gas is assumed to 
carry a very large angular momentum, so that a disk forms with an outer radius at more 
than 10 5 Schwarzschild radii. To suppress the overall emission, which now includes the 
additional dissipation of this large angular momentum, it is also assumed that the electron 
temperature is much lower than that of the protons (T e <C T p ). In fact, T e < 10 10 K. Since 
the electrons do the radiating, the efficiency remains small even though the protons are 
very hot. 



1.2. A Sub-equipartition Magnetic Field 

The idea that Sgr A*'s low radiating efficiency is due to a sub-equipartition magnetic 
field B deserves further attention, especially in view of the fact that the actual value of B 
depends strongly on the mechanism of field line annihilation, which is poorly understood. 
Two processes that have been proposed are (i) the Petschek (1964) mechanism, in which 
dissipation of the sheared magnetic field occurs in the form of shock waves surrounding 
special neutral points in the current sheets and thus, nearly all the dissipated magnetic 
energy is converted into the magnetic energy carried by the emergent shocks; and (ii) van 
Hoven's (1979) tearing mode instability, which relies on resistive diffusion of the magnetic 
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field and is very sensitive to the physical state of the gas. In either case, the magnetic field 
dissipation rate is a strong function of the gas temperature and density, so that assuming a 
fixed ratio of the magnetic field to its equipartition value may not be appropriate. 

Kowalenko & Melia (1999) have used the van Hoven prescription to calculate the 
magnetic field annihilation rate in a cube of ionized gas being compressed at a rate 
commensurate with that expected for free-fall velocity onto the nucleus at the Galactic 
Center. Whereas the rate of increase dB/dt\f in B due to flux conservation depends only 
on the rate f of the gas, the dissipation rate dB/dt\d is a function of the state variables and 
it is therefore not necessarily correlated with f. Although these attempts at developing a 
physical model for magnetic field dissipation in converging flows is still rather simplistic, it 
is apparent from the test simulations that the equipartition assumption is not always a good 
approximation to the actual state of a magnetohydrodynamic flow, and very importantly, 
that the violation of equipartition can vary in degree from large to small radii, in either 
direction. As such, calculations that assume equipartition of the magnetic field with the 
gas throughout the domain of solution may be greatly underestimating the importance of 
the deviations of B from its B eq value since the predicted spectrum relies critically on the 
contribution from magnetic bremsstrahlung. 

1.3. Impact on Sgr A*'s Spectrum 

The first serious attempt at modeling the spectrum of Sgr A* as being due to emission 
by the accreting gas was carried out by Melia (1992, 1994), who assumed a black hole 
mass of ps 10 6 M Q . But this mass is no longer consistent with the now more accurately 
known value of ~ 2.6 x 10 6 M , which accounts for a factor of ~ 7 increase in M . In 
addition, the earlier calculations integrated the cyclo-synchrotron emissivity out to the 
lowest 20 harmonics only, which misses some of the contribution to the flux by the highest 



temperature gas at the smallest radii (Mahadevan, Narayan & Yi 1996). In this paper, we 
recalculate the spectrum produced by a quasi-spherical infall onto Sgr A* using the updated 
mass value, a more accurate handling of the magnetic bremsstrahlung and an empirical fit 
to the magnetic field, motivated by the simulations of magnetic dissipation discussed above. 
In section § 2, we derive the equations governing this spherical infall, which we adopt as a 
simplified version of the real accretion picture. Of course, the real accretion flow will deviate 
from radial at small distances from the black hole, where the gas begins to circularize 
with its advected specific angular momentum. In a fully self-consistent calculation, we 
will use the distributions derived from an actual 3D hydrodynamic simulation as the basis 
for calculating the emissivity. The model parameters and results of our calculations are 
discussed in § 3, and we summarize our conclusions in § 4. 

2. Equations Governing Spherical Accretion 

2.1. The Radial Profiles 

We follow the sequence of derivations in Shapiro (1973), with the primary differences 
being the inclusion of the magnetic field and a radiation pressure term and the fact that we 
restrict out attention to supersonic flows. The equation of mass conservation reduces to 

M = 4vrr V , ( 2 ) 

where M is the mass accretion rate onto the black hole, p is the mass density of the 
accreting gas, — v = u r is the radial component of the fluid 4-velocity (but defined to be 
positive inwards, so that v = — dr/dr for the infalling plasma), and r is the distance from 
the black hole. Equation @ can be recast into the form 
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where a prime denotes d/dr. At Rq, where the numerical integration begins, we assume 
that the gas is non-relativistic and that the gravitational potential is weak. 

The second equation, arising from momentum conservation, is the steady state 
relativistic Euler equation for a spherical geometry: 



,'r? + v 2 -2GM/r\ n , GM 



where the mass of the central black hole is given by M, the total pressure is given by 

P=ff + fl». (5) 

the non-magnetic pressure is given by 

P th = 2nkT + P rad , (6) 

the particle mass-energy density is 

e p = m p c 2 n , (7) 
and the internal energy density of the gas is 

B 2 

e = ankT + 3P rad + —. (8) 

In the fully ionized but non- relativistic limit (i.e., 10 5 < T < 6 x 10 9 K), a = 3. On the 
other hand, in the relativistic electron limit (6 x 10 9 < T < 10 13 K), a = 9/2. We use the 
general expression from Chandrasekhar (1939) that is valid for all T: 

where x = m e c 2 /kT, y = m p c 2 /kT and Kj refers to the i th order modified Bessel function. 
We assume that the gas consists solely of completely ionized hydrogen. Note that since 
we will assume a mostly radial B (see below), the numerator of Equation (f|) does not 
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depend on the magnetic field. That is, there is no large scale current, since v x B = 0. The 
radiation pressure, P ra d, is given by the Rayleigh- Jeans approximation 



where v m is the frequency below which the radiative emission is highly absorbed, so that 
the optical depth, T^°(u m ), from r to infinity is unity. This assumes that the accretion 
flow is relatively unimpeded by radiation pressure, or, in other words, that the accretion 
is sub-Eddington. As mentioned in Melia (1992), with L Sgr A* ~ 1O 5 L <C L E d ~ 10 n L©, 
(Zylka et al. 1995; Shapiro and Teukolsky 1983), this certainly appears to be the case for 
Sgr A*. Numerically, this means that hu m <C kT, which we verify a posteriori. 

The third primary equation follows from energy conservation, which we derive from 
the first law of thermodynamics: 



where r is the proper time in the gas frame. The use of Pth rather than P assumes that 
the compression of the gas is parallel to the magnetic field lines (i.e., radial inflow with a 
mostly radial B). The heating (r) and cooling (A) terms render the flow non-adiabatic (i.e., 
ds/dr 7^ 0). The radiative cooling includes magnetic bremsstrahlung and electron-ion and 
electron-electron thermal bremsstrahlung (see Melia 1994, with an updated prescription 
in Melia & Coker 1999). Local UV heating from nearby stars results in a minimum gas 
temperature of 10 4_5 K (Tamblyn et al. 1996), but the shocked gas at the model's outer 
radius, Rq, is expected to be hotter. Thus, the major heating term (other than the effects 
of compression) is expected to be due to magnetic field reconnection. Specifically, we use 
(see Ipser & Price 1982) 




(10) 




(11) 




(12) 
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Thus, if the magnetic field is flux conserved, for which B(r) oc r~ 2 , then no reconnection 
is taking place and T = 0. However, if the magnetic field remains in approximate kinetic 
equipartition with the gas, meaning that its pressure increases in tandem with the ram 
pressure of the accreting gas (P mag / Pram ~ constant), then B(r) oc r~ 5//4 , and r ^ 0. Since 
T(r), v (r), and n(r) will turn out not to be perfect power-laws, these statements are only 
approximate. 

To find the velocity gradient, v', we place the above definitions for P th and P ra d into 
Equation (f|) to get 

( c 2 + v 2 -2GM/r \ f 8k7c 3 2 , \ GM 

vv =-{ P + e p + e J [ 2k ^ nT + nT ^ + ^ Tiy rn + T3iy m u m }J - — ■ (13) 

Thus, substituting for n' from Equation (13), we have, 



2 



, _ -Hv [2nk{V - 2T/r) + (8fc7r/9c 3 )(TV^ + T3z/>' m )] - vGM/r 
V ~ v 2 - 2nkTH ' { } 

where, for ease of writing, we have defined the quantity 

* + *-2GM/r 

P + e p + e K J 

Note that in the non-relativistic and small B limit, H = 1/p. For simplicity, we write v' as 
the sum of two terms: 

v' = fT' + g, (16) 



where 



and 



-Hvk(2n + 8'Kui/9c 3 ) 
1 v 2 -2nkTH ' 1 ' 



-HvkT (-An/r + 8?rz/>' m /3c 3 ) - vGM/r 2 



v 2 - 2nkTH 

The form of these expressions is that of the classic wind equations (see, Parker 1960; 
Melia 1988). In the simulations we consider here, the gas is supersonic at Rq and remains 
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supersonic on its inward trajectory (i.e., for r < Rq). We, therefore, avoid the special 
handling required for solutions that cross any sonic points, where the denominator of 
Equation (|14}) vanishes. 

Since d(e p /ra)/dr = 0, Equation (|TT|) reduces to 

7T (-) = - p *T (-) + ' 19 

dr \nj dr \nj n 

Substituting in P t h and e, we get 

d / , m MvlkT B 2 \ ( lrTi MvlkT\-n' T-A 

~ "dF [ akT + iZr + 8^) = {- 2nkT ~ J V + — ■ (20) 

Substituting in Equation ([3D and explicitly taking the derivatives, we then obtain 
T [ ak + „ "! = - — + - 2A;T + — + 



3nc 3 y U r / V 9nc 3 87m y 

A-r gg 8 7 r^fcTz/ m / 
m> 47rn nc 6 

We use Equations ( fH|) and fl2"T|) to determine v' and T', respectively, then solve for v 
and T using an implicit differencing scheme, with 

T j+1 = Tj + (r j+1 - r )T' J+l , (22) 

where T'- +1 is a function of Tj,Vj, Tj, Bj, and v'y Similar relations hold for v. Thus, the whole 
flow is determined given some outer boundary conditions (Rq, (3, f v , and fx] see below), 
the empirically inferred behaviour of B(r), and a prescription for finding u m . Note that 
boundary condition values for T' , v', B', and v' m are also needed; at the outer boundary, we 
assume the flow is optically thin, in adiabatic free-fall with a thermal equipartition field 
and determine the derivatives accordingly. 



-14- 



2.2. The Optical Depth and a Prescription for u m 

Since the gas is not expected to be optically thin at low radio frequencies, the effective 
optical depth as a function of frequency must be determined before the radial profiles can 
be calculated. For Sgr A*, we assume (and check a posteriori) that the energy absorbed at 
any given radius is small compared to the thermal and kinetic energy at that radius; this 
is not likely to be true for objects that are accreting closer to their Eddington luminosity. 
Further, we assume that multiple scatterings are unimportant. 

Following Rybicki & Lightman (1979), we use, for an effective optical depth in a zone 
of observed size dr (= r^+i — rj) at infinity, 



Tj = Cdr^J a abs (a abs + na scat ) } (23) 

where 

dr = dr 1-/3/1 , (24) 

^(1 -2GM/rc 2 ){\- (3 2 ) 

a a bs is the absorption coefficient, n is the electron number density, and o sca t is the electron 
scattering cross section. Note that r J= i = lr s = 2GM/c 2 . The coefficient C in Equation 
( p3|) is a geometric term arising from the fact that in spherical symmetry, the average path 
length of a photon that reaches the observer is somewhat larger than drQ. It is given by 

C = mui{-L, — } , (25) 

Z 1 Umax 

where fi max is defined in Equation (p7|). In Equation fl24|), 

(3 = (26) 
cJl + (v/c) 2 - 2GM/rc 2 



and /i is the average of the minimum of the cosine of the angle between the line of sight and 
the flow (which is here —1) and the maximum, 
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given by Zeldovich & Novikov (1971). This maximum in /x is a consequence of the fact 
that not all of the emitted photons reach the observer; some are captured by the black 
hole. Note that p max changes sign and becomes negative at radii smaller than 1.5 r s . Also, 
fJ'max —1 at the event horizon since the curvature in the trajectory of photons emitted 
at other angles takes them back into the event horizon. Although r has a more complex 
angular dependence, the use of p, in Equation (E3) is a necessary simplification at this point. 



For a a t, s , we use Kirchoff's Law, 

a abs = 3u/B u , (28) 

where j v is the total emissivity (in ergs cm -3 s _1 Hz -1 steradian -1 ) and B u is the blackbody 
Planck function. For the models discussed here, we use a total emissivity that includes 
cyclo-synchrotron emission (Coker & Melia 1999) and electron-ion and electron-electron 
bremsstrahlung (Melia & Coker 1999). For a scat we use the exact scattering cross section 
(see, e.g., Lang 1980) 



3 

O'scat — -^&T 



1 + x (2x(l + x) , , A 1, . , l + 3x 
1 ln(l + 2x) + — ln(l + 2x) 



(29) 



x 3 \ l + 2x v 7 2x v ' (l + 2a;) 2 

where x = , yhi , /(mc 2 ) and 

7 = max{l, Vl2kT/mc 2 } (30) 
is the RMS value of the thermal velocity at the temperature in the given zone. 

Finally, to find the total optical depth from zone j out to infinity at some observed 
frequency u , related to the emitted frequency v by 

v /(l-2GM/rc 2 )(l-/5 2 ) 

Uo = u T=ft ' (31) 

we use 

k=oo 

C = E (rj/nfn . (32) 
k=j+l 
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The minimum frequency v m that a photon needs to have in order to escape is found by 
iteratively determining the wavelength at which r k is unity with the caveat that u m not be 
less than the plasma frequency 

v p = ej-^— . (33) 
y nm e 

This caveat is required since photons with a frequency less than v v are unable to propogate 
and are thus trapped by the infalling gas. Once the radial profiles (optical depth, density, 
velocity, and temperature) are determined, it is possible to calculate the emission spectrum 
for a given magnetic field profile. 



2.3. Calculation of the Spectrum 

Now we are ready to calculate the predicted observable luminosity L uo at infinity (see 
Shapiro 1973; Ipser & Price 1982; and applied to Sgr A* in Melia 1992, 1994): 

where, if the emitting zone is optically thick (i.e., if Tj > 1), then 

l v = B V [}-e~ ri ] ■ (35) 

Otherwise 

I v = j v e- T 'Ar. (36) 

In these two expressions, B v and j v depend on \i not p,. The use of Tj here is a calculational 
simplification; ideally, it too should be a function of /i. The sum over j is truncated at J, 
for which rj = R . It is assumed that r^°(J) = 0. This ignores the possible absorption by 
Sgr A West of the low frequency (i>q < 10 9 Hz) radiation (Beckert et al. 1996). Sgr A West 
is an HII region surrounding Sgr A*. 



3. Parameters for a Spherical Accretion Model for Sgr A* 



The stellar winds that accrete onto Sgr A* are thought to be hypermagnetosonic, so it 
is sensible to continue using the standard definition of an accretion radius (Eq. [p]]), which 
is defined to be the point at which the gravitational potential energy is equal to the initial 
kinetic energy of the gas, but using an average of the stellar wind velocities weighted by 
the mass loss rate. For the Galactic Center, Ra ~ 10 16 ~ 17 cm. However, this definition of 
Ra does not necessarily give the physical radius at which the capture always occurs. For 
example, wind-wind collisions, combined with radiative cooling, can reduce the effective 
value of v w . For simplicity, we use this standard definition of Ra to set the primary length 
scale of the problem, and put Ro = /rRa, where fn is a parameter of order unity. 

Now, M is likely to be a fraction of the total GC wind, M w , which is estimated to be 
~ 3 x 10~ 3 Mq yr _1 (Geballe et al. 1991). However, more recent work suggests that the 
mass loss rate of at least some of the central massive stars is less than previously thought 
(Hanson et al. 1998), so that M w may be smaller as well. If the GC wind is dominated by 
a source some distance D > Ra away from Sgr A*, then one would expect 

• _ R\M W _ (2GM) 2 M W 
~ D 2 ~ v^D 2 ■ {67) 

For example, if the 1000 km s _1 wind from IRS 13E1 (which is thought to be ~ 0.16 pc 
away from Sgr A*; see Melia & Coker 1999), dominates the flow, then M ~ 10 21 gs _1 . 
Observationally, the mass accretion rate is not very well constrained; a reasonable range 
is ~ 10 20_22 gs _1 . For all of the models presented here, we assume that the mass M of 
the central black hole is 2.6 x 10 6 M Q , based on the latest observations (Eckart & Genzel 
1999; Ghez et al., 1998). The magnetic field B at 1 Ra is thought to be a few milliGauss 
(Yusef-Zadeh et al. 1996) but its large scale average value may be a few times smaller than 
this (Marshall et al. 1995). In practice, we use a value for B such that at Rq the magnetic 
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field is close to thermal equipartition: 

B 2 /8n = P mag = r]P thermal = r]2nkT, (38) 

where 77 is a parameter of order unity. Together, M and 77 effectively determine the 
normalization of the model spectrum. 

^From a combination of Hel line observations (Najarro et al. 1997) and hydrodynamic 
modeling (Coker & Melia 1997), the accreting wind is thought to be supersonic at _R with 
a temperature T ~ 10 6 ~ 7 K and a velocity Vq ~ 10 2 ~ 3 km s _1 . The temperature, T w , and 
velocity, v w , of the individual stellar winds that accrete onto Sgr A* are not necessarily the 
same as T and vq. Wind-wind shocks and subsequent radiative cooling and acceleration 
result in T ^ T w and vq ~ v w . The winds, which originate well outside of 1 Ra, are 
originally hypersonic with a Mach number of 10-30, but by the time they reach 1 Ra, they 
are on average only mildly supersonic with a Mach number of 1-3. The non-zero velocity at 
infinity tends to make Vq > Vff, where v // is the free-fall velocity, while shocks tend to have 
the opposite effect, making vq ~ f///4. Therefore, for simplicity, we assume that vo is given 
by a fraction f v , of order unity, of the free-fall value. With R = /rRa and Equation ([[]), 
we then have 

l 2GMf R 

vo = fvVw = /«y — — • ( 39 ) 

In order to end up with a supersonic accretion solution, it is also necessary to set To to be 
less than the local virial temperature, i.e., 

2Wo < GM_ ^ 2kT mr 
m p 2R m p 

We parameterize the starting temperature according to the prescription T = fxT vir , where 
again fx is of order unity. 
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3.1. Angular Momentum Considerations 

Some accretion models of Sgr A* (e.g., Narayan et al. 1995) assume that v <C i>//, 
requiring a deceleration mechanism fairly close to 1 Ra- Since T and B are too small to 
decelerate the inflowing gas, the most likely reason for vq being so small would probably be 
that the gas contains a greater specific angular momentum than is seen in the hydrodynamic 
simulations (Ruffert & Melia 1994; Coker & Melia 1997). With the wind originating from 
more than a dozen sources distributed fairly isotropically around Sgr A* (Melia & Coker 
1999), the time averaged specific angular momentum accreted by the black hole is likely to 
be small. This is the conclusion to which one is led with the current series of hydrodynamic 
simulations. The circularization radius at which the accreted specific angular momentum is 
equal to its Keplerian value is 

R circ = 2\ 2 r s , (41) 

where r s = 2GM/c 2 is the Schwarzschild radius and A is the accreted specific angular 
momentum in units of cr s . The hydrodynamic simulations suggest that A is highly time 
variable with a value in the range 2 — 50. However, these calculations have had limited 
physical resolution near the black hole and they have therefore sampled the A accreted 
within the region Rboundary r s . As one would expect, tests have shown that smaller values 
of Rboundary result in smaller values of A carried across that threshold. Thus, the average 
values for A given above may be thought of as upper limits so that R C i rc <C Ra- Similarly, 
the viscous transport of angular momentum outward and a possible mass loss via a wind 
off the disk (see, e.g., Xu & Chen 1997) should result in M(R circ ) > M(r s ). Numerical 
calculations as yet do not have the physical resolution to verify this in the case of Sgr A*. 

More importantly, these simulations indicate that the vector direction of the accreted 
specific angular momentum is highly time variable and dependent on conditions outside 
of the accretion radius, with sign changes occurring on a time scale of 1-50 years. Thus, 
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with R C i rc representing the length scale within the nascent disk (though the distribution 
of angular momentum outwards may increase the disk's outer radius somewhat above this 
value), even if such a disk forms, it is not likely to be stable over a time scale of more than a 
few decades, since the gas contains no long term preferred angular momentum axis. Recent 
work (Genzel 1998) suggests that the wind sources may have a small net rotation around 
Sgr A* and, if the resulting additional angular momentum is not canceled in wind-wind 
collisions, then angular momentum could play a more significant role in the spectrum of Sgr 
A* than we assume here. 

3.2. The Magnetic Field 

In this paper, we do not attempt (as Kowalenko & Melia 1999 did) to determine 
the magnetic dissipation rate from first principles, but rather we will use our fits to the 
observed spectrum of Sgr A* to infer empirically what the profile of B should be in order 
to validate this picture. In other words, we have a sense of how B might behave (most 
likely different from simple equipartition) based on the calculations of Kowalenko & Melia 
(1999), and we use that algorithm for B{r) to fit the spectrum. We find that a 3 component 
profile is necessary and sufficient for this purpose. The profile of B affects primarily the 
heating term, T, and the magnetic bremsstrahlung emissivity, jl ync , as discussed above. 
With 3 components, we have 5 parameters describing the magnetic field: 3 power-law 
indices (pi,P2,P3) and the location of 2 breakpoints (ri,r 2 ). Since the magnetic field of the 
ISM is observed to be approximately in thermal equipartition, we fix the first index, p 1 , 
so that B{r) oc r~ 5 / 4 , or P mag = i]P ram . The shape of the final spectrum is sensitive to 
the remaining 4 parameters (ri, p 2 , r 2 , and p 3 ) but we find that one general configuration 
works best. The first breakpoint r\ occurs at a few hundred r s , with B(r) — > constant 
(i.e., P2 — 0). Interestingly, the second breakpoint r 2 is at ~ 3r s , the radius of the last 
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stable orbit for a non-rotating black hole. At that point, B(r) oc r 3 6 , which suggests the 
presence of a magnetic dynamo at small radii. 

These three segments of B(r) (an equipartition region, a region of constant field 
intensity, and a dynamo region) are all necessary to fit the general flat spectral shape 
of Sgr A*. The equipartition region results in a sufficiently large field and temperature 
to reproduce the observed flux at low frequencies (i.e., 10 9 ^ v <, 10 10 Hz) with r large 
enough so that the emission is not self-absorbed. The region of constant magnetic field is 
required to match the flat spectral index of ~ 0.3 (Lo 1986) at moderate radio frequencies 
(10 10 v ;$ 10 11 Hz). At ~ 2 x 10 11 Hz there is an observed sub-millimeter excess, 
presumably due to an ultra-compact component a few r s in size (Falcke et al. 1998); our 
dynamo region produces this excess. There are many arguments (see, e.g., Hawley & Balbus 
1992) for the presence of a magnetic dynamo close to a black hole so such a component is 
not implausible. 

The calculation of the magnetic bremsstrahlung contribution to j u assumes that /i, the 
cosine of the angle between the flow and the line of sight, is the same as the cosine of the 
angle between B and the line of sight. For a spherical flow, in the absence of magnetic 
reconnection (which should be minimal for r ^> Ra anyway), B r oc r~ 2 , while Bg oc r _1 . 
Thus, even if the magnetic field at infinity is perfectly tangled, it is reasonable to assume 
that at i? , B r 3> B e . However, because the magnetic field is divergence-less, the stretched 
out field lines must close somewhere. Presumably this occurs quite close to the black hole. 
In a fully self-consistent calculation, we will use the distributions derived from an actual 3D 
hydrodynamic simulation as the basis for calculating the emissivity. 
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4. Results 

We present the resulting spectrum of our best fit model in Figure |l|(a). For comparison 
with previous calculations (e.g., Melia 1994), we also plot results for a model which assumes 
only gravity acts to accelerate the accreting gas (i.e., H = 0). Since the flow is only mildly 
supersonic at moderate radii, the assumption of free-fall is physically invalid, but such 
models do allow hotter gas to exist at larger radii and, thus, have more emission at lower 
frequencies than the full solution models. Nonetheless, both model spectra are consistent 
with the observations over 16 decades of frequency. The parameters used in the two fits are 
given in Table |l|. Note that one clear prediction of these models is that the source of the 
majority of the high energy 7-ray emission is not Sgr A* or its associated accretion flow, 
unless some exotic particle acceleration mechanism is operating within the inflow. This is 
consistent with recent models (Melia, Yusef-Zadeh, & Fatuzzo 1998), which suggest the GC 
high energy emission is due to Sgr A East, a super- nova-like remnant near or possibly even 
enveloping Sgr A*. The details of the predicted emission in the IR region are fairly sensitive 
to the details of radiative transport, since the gas circularizes at small radii, possibly 
producing a small disk which emits primarily in the IR. Nonetheless, the next generation 
of infrared observations should be able to detect the emission from the inner ~ 10r<, of the 
black hole. 

Some of the observations used in Figure |T](a) are averaged over a span of years while 
others are contemporaneous. Since the flux of Sgr A* can vary with time by more than 
25% (depending on frequency), we do not necessarily wish to reproduce the entire suite of 
observations but rather the general trend. The resolution of the observations also varies 
greatly but most of the values reported in the table have backgrounds subtracted out and 
in principle refer to the flux of Sgr A* only. Those observations which have a potential 
for source confusion are plotted as upper limits. For example, at v ^ 8.6 x 10 11 Hz, the 
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Table 1: Parameters for the Sgr A* Spectral Fits 



Parameter 


Full Solution 


Free-fall 


Moyr-^gsec- 1 ) 


8xl0 21 


5xl0 21 


Woo (km s _1 ) 


500 


1000 


/r 


1.0 


1.0 


St 


0.5 


1.0 


fv 


0.5 


1.0 


V 


0.05 


0.15 


Pi 


2.5 


2.5 


P2 


0.0 


0.0 


P3 


6.0 


12.0 


n(r s ) 


750 


2000 


r 2 (r s ) 


4 


3 



Model parameters for two fits to the spectrum of Sgr A*, assuming 
spherical accretion. The full solution is for a calculation using the equations 
described in the text while the free-fall model assumes if = so that only 
gravity accelerates the infalling gas (see Eq. ||14|| ). 




Fig. I. — The solid curves are for a model which assumes pure free-fall (H = 0) while the 
dotted curves are for a model which uses the complete equations (see text for details), (a) 
The observed spectrum of Sgr A* along with the predicted results for the best model fits. 
The observational data are taken from various sources (see Narayan et. al. 1998 for a recent 
compilation), (b) The observed size of Sgr A* along with the predicted size, defined as the 
radius at which the optical depth, r™{r) equals 2/3. 
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background emission due to, e.g., dust at a temperature of 90K, is comparable to emission 
from Sgr A* itself (Zylka et al. 1992). This is particularly true in the infrared, where Sgr 
A* has yet to be definitely identified and has a flux considerably less than that of the 
surrounding gas and stars, and in the X-rays and 7-rays, where the poor spatial resolution 
is likely to result in multiple sources being in the field of view. Similarly, observations at 
v 1 GHz are unreliable as limits for the luminosity of Sgr A* due to source confusion, 
scattering, and self-absorption (Davies et al. 1976). In converting the observed fluxes to the 
luminosities plotted in Figure |l|(a), it is assumed that Sgr A* is 8.5 kpc away. 

The predicted size of Sgr A* is shown in Figure [1](b). Also shown are the present 
observational limits and measurements (see Lo et al. 1998 for a summary of recent work). 
At present there are two observations which address the intrinsic diameter of the minor 
axis of Sgr A*. However, these observations are difficult and, for example, may have had 
difficulties with calibration (see, e.g., Krichbaum et al. 1998). The model satisfies the 
upper limit at 3.5mm as well as the lower limit at 0.8mm, but it predicts a source that is 
somewhat more compact than the other observations suggest. But while the observations 
tend to fit Gaussian FWHM to the data, the predicted diameter is defined as the last 
scattering surface (i.e., r[r = 2/3]). For a more detailed and accurate size comparison, one 
needs a more sophisticated treatment, using, for example, the theory of wave propagation 
in an extended, irregular medium (Melia, Jokipii & Narayanan 1992). In addition, at small 
radii, the presence of angular momentum and asymmetries in the flow will result in a 
distended non-circular source shape. For example, at 7mm the axial ratio is observed to be 
less than 0.3 (Lo et al. 1998). Thus, the model size results presented here should be taken 
as approximate lower limits. 

In Figure |2| we show the radial profiles for the temperature and magnetic field for 
the two fits shown in Figure [l]. The gas is always sub-virial (T vir ~ 10 12 K [r/r s ] -1 ) with 



-26- 



a peak temperature of ~ 10 11 K. This peak temperature is within the present observed 
lower (1.3 x 10 10 K; Lo et al. 1998) and upper (5 x 10 n K; Gwinn et al. 1991) brightness 
temperature limits. Only within ~ 10 2 ~ 3 r s of the black hole, where the gas becomes 
relativistic, is there significant magnetic bremsstrahlung emission. It is not a coincidence 
that r\ is ~ 10 3 r s . The accreting gas must stay hot enough at such radii to produce the flat 
spectrum; the constant magnetic field and subsequent heating via reconnection very nearly 
achieve this. However, this cannot continue all the way down to the event horizon or else 
the resulting magnetic field is insufficient to produce the observed sub-millimeter excess 
and the gas gets too hot to continue accreting and becomes pressure supported. Thus, 
we require a sharp increase in the magnetic field strength within ~ 10r s of the black hole 
with a resulting peak magnetic field of ~ 100 Gauss. We postulate that this may occur 
when the frozen-in magnetic field lines, which have been stretched out during the accretion 
process, finally close or are twisted by sheared gas motions due to a transition resulting 
from residual angular momentum in the flow; however, if the viscosity is large, the bulk gas 
flow still remains primarily radial and supersonic (Narayan et al. 1997). As described by 
Hawley and Balbus (1992), instabilities in the flow, possibly tied to differential rotation 
once the gas circularizes, may result in the generation of strong poloidal and toroidal fields, 
driving the magnetic field strength up to thermal equipartition, as required in our model. 
Except in the transrelativistic region (~ 10 3 r s ), the gas density and velocity profiles are 
close to that from free-fall so that n~ 10 10 (r s /r) 3 ^ 2 cm~ 3 . 

5. Conclusions 

By incorporating an enhanced treatment of the magnetic bremsstrahlung emissivity 
and solving the accretion flow equations explicitly, we have improved the accreting black 
hole model for Sgr A* with resulting spectra that are consistent with observations over more 
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Fig. 2. — Plots of temperature (a) and magnetic field (b) versus radius for the two best fit 
models. The solid curves are for a model which assumes pure free-fall (H = 0) while the 
dotted curves are for a model which uses the complete equations (see text for details). 
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than 16 decades of frequency. The mass accretion rate of ~ lCr 4 M yr _1 , as determined 
from the spherical accretion model, is consistent with the rate expected on the basis of 
observations and hydrodynamical arguments (see, e.g, Coker & Melia 1997). However, other 
models, such as an "advection dominated accretion flow" (ADAF) with outflow (Blandford 
& Begelman 1999), require an accretion rate that is 2 or more orders of magnitude smaller 
than this. 

There are some difficulties with the model. The full solution, using the Euler equation, 
suggests that the flow may be transonic and that it may therefore shock at smaller 
radii. A shock could produce high energy particles, which would result in significant 
magnetic bremsstrahlung emission at larger radii (Markoff, Melia, & Sarcevic 1997), thereby 
increasing the emission at lower frequencies. We hope to address this difficulty in future 
work. It is interesting that the break seen in the low energy 7-ray observations coincides 
with the high energy turnover in the model, suggesting the presence of another, hotter 
source within the field of view of EGRET. In addition, we have effectively ignored the role 
played by angular momentum, the presence of which would introduce viscous processes 
that are even more poorly understood than magnetic reconnection. Also, the handling of 
radiative transport has so far been very simplistic. Inclusion of relativistic 3D radiative 
transfer, including ray bending and inverse Compton, could potentially alter the resulting 
spectra significantly. 

The model suggests that if the accretion rate onto Sgr A* is large (~ lCr 4 M yr _1 ) 
then the observed X-ray emission from the GC is due to the extended X-ray emission 
from this accretion process. However, the observations have a resolution on the order of 
arcminutes while the 'extended' emission region from the model spans a few arcseconds. 
Due to calibration difficulties and an uncertain column depth between here and the GC, 
the observed X-ray limits are uncertain to within a factor of a few (Narayan et al. 1998) so 
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the model is still marginally consistent with the X-ray observations. The next generation of 
observations (such as with Chandra) will most likely settle this issue. 

The picture with which we have worked here under-predicts somewhat the low 
frequency radio emission. Within the context of the model, unreasonable temperature and 
magnetic field profiles are required to fit the observations and even then the fit at other 
frequencies degrades. This difficulty is due to two things. First, the low frequency emission 
from close to the black hole is trapped so it does not contribute to the low frequency 
spectrum. Second, gas with sufficient magnetization and temperature to produce the 
observed magnetic bremsstrahlung emission at large radii is gravitationally unbound and 
could not self-consistently accrete. An improvement to the model that addresses this low-z/ 
deficiency in the spectrum may come with more realistic 3D simulations, which will be 
reported elsewhere. 

Finally, although we have briefly touched on the likely importance of the sub-millimeter 
excess to our understanding of the environment just outside the event horizon, this 
clearly is an issue that warrants further detailed theoretical work. Taking into account 
scatter-broadening of the image in the interstellar medium, and the finite achievable 
telescope resolution, the ~ 10r s "shadow" of the Galactic Center black hole (see, e.g., 
Falcke, Melia & Agol 1999) should be well observable with very long-baseline interferometry 
at sub-millimeter wavelengths. In our picture, the dynamo effect that leads to the 
intensification of the magnetic field that accounts for this spectral excess may be due to the 
circularization of the infalling gas when it approaches the circularization radius 2A 2 r s (see 
§ 3.1). We are in the process of examining the behavior of the magnetized plasma when 
shearing motions become important in this region, and we will report on the results of this 
investigation in the near future. 
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